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I Three different hybrid Vlasov-fiuid systems are derived by applying reduction by symmetry 

^ ' to Hamilton's variational principle. In particular, the discussion focuses on the Euler-Poincare 

I formulation of three major hybrid MHD models, which are compared in the same framework. 

' These are the current-coupling scheme and two different variants of the pressure-coupling scheme. 

I The Kelvin-Nocther theorem is presented explicitly for each scheme, together with the Poincare 
invariants for its hot particle trajectories. Extensions of Ertcl's relation for the potential vorticity 

' and for its gradient are also found in each case, as well as new expressions of cross hclicity invariants. 
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1 Introduction 

1.1 Hybrid Vlasov-fluid models in plasma physics 

Hybrid Vlasov-fluid plasma models contain elements of both continuum fluids and phase-space proba- 
bility density. The latter obeys a Vlasov kinetic equation, which in turn is coupled to the momentum 
equation of the background magnetized fluid. 

These hybrid kinetic-fluid models arise in several circumstances in modern plasma physics research, 
ranging from fusion research [3G] to astrophysical plasmas [46]. These circumstances involve the 
coexistence of a cold fluid component with an ensemble of energetic particles that require a kinetic 
description. In meeting the challenges presented in such situations, the plasma simulation community 
has begun developing multiscale fluid plasma models that allow hybrid descriptions of the two types 
of flows. This hybrid approach successfully couples ordinary fluid models for the cold fluid component 
to appropriate kinetic equations that govern the phase-space distribution of the energetic particle 
species. 

One research direction of relevance in applications is the development of hybrid schemes for mag- 
netohydrodynamics (MHD) [36]. This development has split into two promising approaches: the 
current-coupling scheme [36, 1] and the pressure-coupling scheme [36, 11, 8, 43]. These schemes differ 
in how the fluid equation is coupled to the kinetic equation for the hot particles. 

Recently, the Hamiltonian formulations of a variety of hybrid Vlasov-fluid plasma models were 
developed that led to new theories of either current-coupled or pressure-coupled hybrid Vlasov-MHD 
models [44] . These Hamiltonian formulations cast considerable light on the energetics of hybrid Vlasov- 
fluid plasmas and the relations between their MHD approximations. In particular, the current-coupling 
scheme has been shown to possess a well defined Hamiltonian structure. However, the pressure- 
coupling schemes were found to require additional fluid transport terms in their accompanying kinetic 
equations to preserve their corresponding Hamiltonian structures, which otherwise would have been 
lost. 

While the Hamiltonian picture of plasma models provides a powerful tool for energy-conserving 
properties and may also open the way to relevant stability considerations, the question has remained 
open whether these models allow a Lagrangian variational formulation by Hamilton's principle. La- 
grangian formulations of plasma physics have been successful in several contexts, since they can be 
easily approached by several approximation methods such as averaging or asymptotic expansions. 
Well known results of this approach are contained in the Lagrangian wave theory by Dewar [40] and 
the celebrated guiding center motion by Littlejohn [26]. The key feature of Lagrangian variational 
formulations is that, in contrast to their Hamiltonian counterparts, the application of essentially any 
approximation scheme preserves the intrinsic geometrical properties of the resulting dynamics that 
emerge from its variational structure in the presence of symmetry. In particular, the derivations on 
the Lagrangian side provide a framework in which straight-forward application of asymptotic expan- 
sions would still preserve the fundamental circulation and Lie-Poisson properties of these theories. 
In contrast, asymptotic expansions of the Hamiltonian formulations, for example, require exceptional 
care in preserving the Jacobi identity, while applications of asymptotic expansions directly to the 
equations of motion typically pay no heed to these geometric properties. An illustrative example of 
this phenomenon for fluids was given in Camassa, Holm and Levermore [()] in deriving the "Lake 
equations" and "Great Lake equations".^ 

The above discussion provides a natural motivation for this paper, whose aim is to present the com- 
plementary derivations of the hybrid fluid models that were introduced in [44] from the Hamiltonian 

^ It was a very good moment in all our careers when Roberto Camassa, Dave Livermore and DDH realized the efficacy 
of applying asymptotics to Hamilton's principle, while working on the board together one afternoon in Los Alamos. 
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side, by recovering them on the Lagrangian, or Hamilton's-principle side. As in the Hamiltonian Lie- 
Poisson formulation, the Euler-Poincare approach presented in this paper incorporates the geometric 
properties that follow from the relabeling symmetry shared by all continuum systems. In addition, 
the Euler-Poincare variational framework provides a systematic framework for the derivation of other 
approximate models, which also inherit these geometric properties from the variational structure. 

The new information we gain in this paper in each case is the natural formulation of a Kelvin- 
Noether circulation theorem and a corresponding Ertel theorem for the potential vorticity. While the 
comparisons of the Lie-Poisson Hamiltonian properties in paper [44] afford insight into the energetics 
of these theories, the comparisons of their complementary derivations on the Lagrangian side provide 
distinctions in their circulation laws, and in their Ertel relations for evolution of the potential vorticity 
and its gradient. Thus, the Euler-Poincare approach affords additional insights into the diagnostics 
of basic MHD processes in the presence of a hot particle Vlasov component. 

Maxwell- Vlasov plasmas have been treated earlier using the Euler-Poincare variational approach 
based on applying symmetry reduction to Hamilton's principle [7]. The present work starts with 
the Low Lagrangian [27] and systematically develops a series of approximate Lagrangians for use in 
symmetry-reduced Hamilton's principles for re-deriving the hybrid Vlasov-MHD fluids in [II]. These 
approximate Lagrangians are shown to admit a variety of symmetry reductions that produce variants 
of Kelvin's circulation law for each theory, together with new expressions for the dynamics of their 
cross helicities, some of which are found to remain invariant. 

The resulting Euler-Poincare equations recover the equations in [44] and illuminate the differences 
in the interplay between the geometric structure and circulation mechanisms of both the current- 
coupled and pressure-coupled hybrid Vlasov-MHD models. The dynamics of the potential vorticity 
and its gradient are also explained through appropriate generalizations of Ertel's theorem to the hybrid 
MHD case. These generalizations arise as a direct consequence of the vorticity dynamics produced by 
the Euler-Poincare equations of motion. 

1.2 Plan of the paper and its main results 

The main content of the paper is, as follows. 

1. The remainder of this Introduction reviews the Euler-Poincare construction of the Maxwell- 
Vlasov equations. Section 2 then extends this system to account for the presence of several cold 
fluid components (kinetic-multifluid system). The reduction process and the resulting circulation 
laws are presented explicitly, including the Poincare invariant relations that are now obtained 
from Noether's theorem. 

2. Section 3 considers the MHD limit of the kinetic-multifluid system, thereby formulating the 
Euler-Poincare equations for the hybrid current-coupling MHD scheme. After presenting the 
Kelvin-Noether theorem, Ertel's relation for the potential vorticity is presented, thereby extend- 
ing Ertel's theorem for MHD [15] to a hybrid model. Also, it is shown that the usual expression 
of cross helicity is conserved by this hybrid model. 

3. Section 4 presents the Euler-Poincare formulation of the first pressure-coupling hybrid MHD 
scheme. In this setting, the cold component drives the whole dynamics, so its velocity adds to the 
mean velocity of the hot particles. This property appears geometrically in the semidirect-product 
Lie group structure that generates the Euler-Poincare construction. The Kelvin circulation and 
Ertel potential vorticity relation are derived explicitly, together with a new expression of the 
cross helicity invariant. 



D.D. Holm & C. Tronci 



Hybrid plasma models 4 



4. Section 5 focuses on the second pressure-coupling hybrid MHD scheme. In this context, the 
assumption of a rarefied hot component allows one to neglect the kinetic energy contribution 
of the corresponding mean flow. Then a decomposition becomes necessary to separate the hot 
particle velocity from its mean flow. This decomposition produces a nested semidirect-product 
Lie group structure that fits into the Euler-Poincare construction. Explicit expressions for Kelvin 
circulation, Ertel's theorem and a new cross helicity invariant again result. 

5. Finally Section 6 summarizes our main conclusions and discusses the outlook for future research 
along the present directions. 



1.3 Euler-Poincare formulation of the Maxwell- Vlasov system 

The variational structure of the Maxwell- Vlasov system has been investigated in many different ways, 
starting from the pioneering work by Low [27]. Since then, several variational formulations of this 
system were presented [37, 38, 39, 47], which are mainly based on Eulerian variables. The Low 
Lagrangian, however, involves a mixture of Eulerian and Lagrangian variables. The first variational 
formulation in terms of purely Lagrangian variables appeared in [7], in which the Low Lagrangian was 
modified by the insertion of an extra term. This extra term ties the Lagrangian particle velocity to 
its corresponding Eulerian coordinate, i.e. x(xo,vo) = v(xo,vo). The variational principle was then 
cast into Euler-Poincare form [21], by using a reduction process that takes advantage of the relevant 
symmetry properties of the Lagrangian under the Lie group of diffeomorphisms (smooth invertible 
maps of both physical space and phase space). 

Motivated by the recent results [44] on hybrid plasma models, one may ask whether the variational 
methods developed in [7] would also apply to hybrid Vlasov-fluid systems. This paper shows that 
these methods do indeed apply and they provide a systematic framework in which to develop a 
fully Lagrangian formulation of the hybrid models. The resulting theory is again an Euler-Poincare 
formulation, which naturally inherits all the Lie-symmetry properties of both fluid motion and Vlasov 
kinetic dynamics. These symmetry properties then provide the various theorems for circulation and 
cross-helicity that are derived later in the paper. 

The present section introduces the approach that we shall follow throughout the rest of this paper. 
In particular, we shall review the Euler-Poincare variational formulation [7] of the Maxwell- Vlasov 
system 

f + v.|^ + A(E + vxB).§^=0 (1) 

ot ax m av 

eo/^oi^ = VxB-g^o/v/d3v, =_VxE (2) 



dt ■ - " .A^uy w ^ - , 

eoV-E = g /"/dV, V-B = (3) 



where q is the particle charge and m its mass, while eo and /io are respectively the dielectric and 
diamagnetic constants. Moreover, in the standard notation adopted here, /(x, v,t) is the Vlasov dis- 
tribution on X M'^, while E(x, t) and B(x, t) are the electric fleld and the magnetic fiux, respectively. 

Euler-Poincare approach. The Euler-Poincare approach to the Maxwell- Vlasov system is based 
on an action principle of the type 



6 /'%„(V',^,$,$,A,A)dt = 

J to 



where the Lagrangian L is a functional 



Lfg : TDiff(rAf) xTQ^R 
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depending on the parameter /o € Den(TAf ) belonging to the space of distributions on the tangent 
bundle TM with local coordinates (xo,vo). Here the notation is such that M is the particle con- 
figuration space, V S Diff (TM) is an element of the Lie group of diffeomorphisms of TM and Q is 
the space of electromagnetic potentials A), i.e. Q = C°°{M) x ^}^{M) (where Q^{M) denotes the 
space of differential 1-forms on M). At this stage, the variational principle produces Euler-Lagrange 
equations on Diff (TM) x Q. The explicit form of the Lagrangian reads as [7] 

^h = \^ j /o( |x(xo, vo)|^ + jx(xo, Vo) - v(xo, vo)|^) d^xod^vo 

-Q j /o(^(x(xo,vo)) -i(xo,vo) • A(x(xo, vq))) d^xo dVo 

+ — /iV^ + atAl^dV-— /iVx ApdV, (4) 

where the potentials determine the electric field as E = — V$ — dtA. and the magnetic flux as B = 
V X A. Also, the map ip in the action principle above denotes (x(xo, vq), v(xo, vq)) := V'(xo, vq) and 
the term 

-m / /o |x(xo, Vo) - v(xo, vo)| d xq d vq 

allows v(xo,vo) to be varied independently and enforces x = v. Dropping the above term returns 
precisely the Low Lagrangian [27]. 

At this point, the invariance property of the Lagrangian (4) is such that 

V', i-. A, A) = L^^„^-i o ci>, 6, A, A) =: /(X, ^, 6, A, A, /) 
where we have defined 

X := V; o G X{TM) , / := /o o € Den(rM) (5) 

and X(TM) denotes the Lie algebra of vector fields on TM. Notice that the dependance on the 
identity element if) o has been omitted in the reduced Lagrangian. In this setting, the reduced 
Euler-Poincare Lagrangian I : X{TM) x Den(rM) x TQ — ?• M produces the equations [7] 

9?5X + ^^^=^^(-'-)57' W + ^^^=°' 
d 61 61 _ d 61 61 _ 

di6^~6^~^' di6A~6A~^' ^' 

Here, the symbol £x denotes the Lie derivative along the phase-space vector field X € X(TM), whose 
components are given by 

X(x, v) = (u(x, v), a(x, v)) 

in which (x, v) S TM are the Eulerian position- velocity coordinates and Af = R^, so X G X(M^). The 
symmetry-reduced version of the Lagrangian (4) is 

/(-m|u| -|--m|u — vj —(/<!> -|- (7U • A j dxdv 



+ ? /iV^ + SjAl^d^x- /|V X Apd^x. (8) 
2 J 2/xo J 



The Maxwell- Vlasov equations are obtained upon applying the variations in the above Lagrangian 
and substituting them into the Euler-Poincare equations (6)-(7), as shown in [7]. 
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Outlook. The remainder of the paper apphes the Euler-Poincare approach to the case of hybrid 
Vlasov-fiuid models that commonly arise in plasma physics research. After studying a general Vlasov- 
multifluid system for the interaction of several fluid plasma components with a hot particle species, the 
paper focuses on comparing the Euler-Poincare structures of current-coupling and pressure-coupling 
hybrid MHD schemes. In the latter case, the geometry of the system provides an interesting example 
of how the Vlasov distribution function may be transported by the background fluid, through diffeo- 
morphisms (smooth invertible maps) acting by tangent lifts. Our considerations here are restricted 
to barotropic fluid flows. 

Prom the strictly mathematical point of view, the case of ideal adiabatic flows that transport the 
specific entropy may be obtained by a straightforward generalization. However, from the physical 
viewpoint the role of heat exchange and the effects of an additional advected quantity should lead to 
other interesting effects that we intend to discuss elsewhere. In particular, adiabatic flow effects may 
be especially interesting for hybrid fluid drift-kinetic models, which may be treated in a Lagrangian 
setting, perhaps by using an approach similar to that for oscillation-center theory, as in [10]. For 
example, this problem might benefit from an exploration of adiabatic invariants that arise from av- 
eraging the Hamilton's principle. This is available for Lagrangian theories, but not for Hamiltonian 
theories, which instead would use Lie series methods. A proof of the equivalence of these theories 
would also be interesting. This approach follows ideas that go back to Dewar [1(1], but now have 
been further illuminated by the advent of the Euler-Poincare approach to reduction by symmetry for 
continuum descriptions on the Lagrangian side [21]. 

An early step in this direction was already made by Holm, Kupershmidt and Levermore in [20], 
who studied Poisson maps in the Eulerian and Lagrangian descriptions of continuum mechanics. Many 
of the concepts from that work, particularly momentum maps from canonical phase spaces to the duals 
of Lie algebras, remain just as important in the present work as they were then, but here they are 
applied on the Lagrangian, or Hamilton's principle side for hybrid Vlasov-fluid systems. 



2 Vlasov-multifluid system 

This section presents the Euler-Poincare formulation of a system composed of several fluid plasma 
species, each denoted by its label s = 1, . . . ,N, with an energetic Vlasov component. The Hamiltonian 
formulation of this system has been presented in [44] and its equations of motion are expressed as 

dU 

+ • V) ?7, = asPs (E + UsXB)- Vp, (9) 

ot 

^ + V • (psUs) = (10) 
df df , , df 

^ + v-^ + a/, E + vxB --^ = 11 
ot ax av 



Po^o~Q^ = V X B - ^0 X a-sPsUs - p-oqh rv f d^v (12) 

s 

- = -VxE (13) 
eoV ■B = ^asPs + qh [ fd\, V-B = (14) 

s 

In these equations for the Vlasov-multifluid system, = Qs/mg is the charge-to-mass ratio of the 
fluid species s, while ps and Us are its mass density and velocity, respectively, and p<j is the scalar 
partial pressure of species s. In the above system, the index h denotes the hot particle component, 
while each fluid species s is governed by its own momentum and mass-transport equations. In order 
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to avoid proliferation of indexes, the notation V is relegated to denote only spatial gradients applied 
to quantities on physical space. Then, gradients on phase space coordinates are denoted by V(x,v); 
while partial differentiation of phase-space quantities will be denoted by (9x or dv 

For the case that the fluid component is absent, the Euler-Poincare formulation of the resulting 
Maxwell- Vlasov system was presented in [7]. On the other hand, an Euler-Poincare formulation of 
charged fluids was given in [21]. A combination of these two approaches yields the Euler-Poincare 
formulation of the kinetic-multifluid system. Indeed, we shall show that the equations (6)-(7) can 
be suitably generalized to apply for the Vlasov-multifluid system. As we shall prove below, this 
generalization results from the following choice of Lagrangian: 



/({[/,},{p,},X,/,$,6,A,A) JpslUsfd^^-Y, j Ps{U{ps) + - asU s ■ A) d^ 

+ mh j f(i^ + ^ |u - v|2 - ah^ + ahu ■ A^ d^xdV 



3x 



+ — / |V$ + atArd'^x /iVxArd^x, (15) 

2 J 2po J 

where U{ps) denotes the total internal fluid energy, related to partial pressure of species s by = 
p1U'{ps)-, while the notation {Us} means that that one considers all species s = 1 . . . N (and anal- 
ogously for {ps})- In the special case of a single fluid species, s = 1 and the above Lagrangian is 
defined as a functional 

I : (X(m3) e X (C~(m3)* X C°°(m6)*) x TQ{W^) M 

where X(M") denotes the Lie algebra of vector fields in M", the asterisk denotes the distributional 
dual space and the tangent space TQ(M'^) is constructed on the space Q(M'^) of electromagnetic 4- 
potentials (<I>, A). In this setting, the advected fluid quantity is the mass density p(x) € Den(]R'^), 
while the advected phase-space quantity is the Vlasov distribution /(x, v) G Den(M®). 

At this point, in order to use the above Lagrangian, equations (6)-(7) must be adapted to the 
present case by extending them to account for the presence of the fluid components. In following the 
treatment in [21, 7], one writes the Euler-Poincare theorem in the following general form. 

Theorem 1 (Euler-Poincare kinetic-multifluid system) The kinetic-multifluid system (9) -(14) 
arises from the Euler-Poincare variational principle 

6 r/({C/s},{Ps},X,/,c^,ci>,A,A)dt = (16) 

J to 

with the Lagrangian given in (15) and the variations 

6Us = dtWs-£u.Ws, 6X = dtZ-£xZ, df = -£zf, 6ps = -£wsPs 

where € X(]R'^), Z G X(M^), 5^ and 6A all vanish at the endpoints. This variational principle is 
equivalent to the Euler-Poincare equations 

d SI n SI ^ SI 

d 61 61 61 

^ + £usPs=0, ^ + £^f = o (19) 

I6l_6l_ d_6l_6l_ 

dt§^ 5$"^' dt§A 6A~^- ^^^^ 
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Proof. The equivalence between the variational principle (16) and the Euler-Poincare equations 
(17)-(20) follows easily upon repeating the same steps as in [21, 7]. In order to derive the equations 
(9)-(14), one simply computes the functional derivatives of the Lagrangian (15). In particular, upon 
writing X = (u, a), for the Vlasov kinetic part one has 

61 „ , . , 61 51 mh , ,9 ruh , ,9 . , 

— = m/,/ 2u-v + aftA , ^=0, — = + - + qtu ■ A - qt^ . 

6u oa 0/2 2 

Then, dividing equation (18) by / yields 

d fl 6l\ ^ fl 6l\ ^ 51 , 

diyi^)^^'' J = ^(^'^^ ^7 ■ ^^^^ 

Next, projecting onto the second component yields 



= dvU • dv77 = u(x, v) - V 

/ ou 6} 



so that 



-^(x, v) = (v,a(x, v)) , — = j [mhv + qhA) , — = — -|v| + qh^r ■ A - qh<P . 

ou dj 2 

Upon denoting the particle momentum as p(x, v) := m/^v + (^/^ A(x) and projecting equation (21) 
onto its first component, we obtain 

dtp+ {v ■ + a - 9v)p + 9xV - p = 9x(^|vp + qhY ■ A - qh^^ , 

where (x, v) are independent coordinates. Standard vector identities then produce the Lorentz force, 

/ dA\ 
a(x, v) = - a/j ( V$ + — j + a/j V X (V X A) . 

Therefore, the second equation of (19) gives the Vlasov kinetic equation in the form 



df ^ df 
at ox 



( dA\ 

f V^> + — j - V X (V X A) 



5v 



The fluid equations follow easily by inserting the appropriate variational derivatives into the Euler- 
Poincare equation (17) and the flrst of (19). For example, one computes 

- PsUs + asPsA, — = ^- --as^ + asU-A. 



5Us 5ps dp, 

Next, upon dividing equation (17) by p, one obtains 

dU ( dA \ 

" + £[7^ + a J — + £[7^ A = -Vp, - a,V$ + a,V(J7, • A) 



where = p1U'{ps) is the scalar partial pressure. Finally, equation (9) arises from the explicit form 
of the Lie derivative operation, by using standard vector identities. Analogous arguments also hold 
for the equations of the electromagnetic potentials. ■ 
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Remark 1 (Euler-Poincare reduction) The above theorem follows by an Euler-Poincare reduc- 
tion process applied to the following unreduced Lagrangian 

^fQ,{pos}i{Vs},{Vs},i^,tp, ^',i,A, A) =^J2 I Po^(^o) \Tis{ao)\'^d^ao 

s 

- ^ / pos{aQ){U{pQs) + as^{i]s{ao)) - asi]s{ao) ■ A(r/s(ao))) d^ao 

s 

+ 2 "T-Zi / /o(xo, vo) (^|q(xo, vo)| + |g(xo, vo) - ^^(xo, vo)| j d xq d vo 
j /o(xo, vo)(^$(q(xo, Vo)) + q(xo,vo) • A(q(xo, vo))) d^xod^vg 



+ — / iV^ + StAPd^r- — / iVxAPd^r, (22) 

wzi/i i/ie notation (q(xo, vq), ■u(xo, vq)) = ■(/'(xojVo) /or a group action ip : TM^ — )• TM^. Indeed, the 
invariance property 

^h,{pos}i{ris}, {^s}, ^, ^, ^, A, A) = ^/oo^-i,{po.o„7^}({^^ ° ^ ° ^"^ ^' ^' ^) 
yields the Euler-Poincare Lagrangian 

l{{Us}. {p.}, X, /, CD, 6, A, A) := i/„oV^-i,{p„,or,-}({^^ « ^7^}, o V'"', ^, A, A) 
wiit/i i/ie notation 

Us = 'ns°v7^, X = i/'0'(/;-\ /Os=posO?y7i and / = /oO'0~i. 

T/iis argument follows easily from the treatment in [21, 7]. Notice that for the case of a single species 
■s = 1, the unreduced Lagrangian is of the form 

Lf^^p^^ : rDiff(IR3) X rDiff(r]R3) xTQ^R, 

which emphasizes the Lie group structure that underlies the Vlasov-multifluid system (9) -(14). 

Remark 2 (Kelvin-Noether theorem for the Vlasov multifluid system) It is easy to verify 
that equations (17)-(18) produce the following circulation conservation laws 

/ (?7s(x,t) + a,A(x,t)) •dx = 0, / (v + a?, A(x, f)) • dx = . 

•i* J -ft (Us) cL* -'Ct(X) 

In the first relation the curve 'jt moves with the fluid flow, while in the second relation the curve 
Ct moves with the phase-space vector field X(x, v) = (v, a(x, v)) . More explicitly, one can write 
7t = Vsit) o 70; for a fixed loop 70, and analogously = V'(^) ° Co- This difference emphasizes the role 
of the Poincare invariant associated with the hot particle motion; see [7]. 

The next sections will consider the Euler-Poincare formulation of hybrid Vlasov-MHD models. In 
particular, the discussion will focus on two main types of hybrid systems: the current-coupling and 
pressure-coupling schemes. 

Remark 3 (The Legendre transform) Notice that the Euler-Poincare Lagrangian (15) is degen- 
erate since 61 /5a = 0. This degeneracy is related to a redundancy in the Euler-Poincare construction, 
which carries all the information about particle paths that are already encoded in the Vlasov equation. 
As explained in [7], this degeneracy presents some problems when one wants to perform a Legendre 
transform to obtain the corresponding Hamiltonian description. However, these problems may be over- 
come by a standard use of Dirac constraints. Similar arguments to those in [7] also hold for the hybrid 
models treated in this paper. 
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3 Current-coupling hybrid MHD scheme 

3.1 Formulation of the model 

In common physical situations, one is interested in single-fluid models. In the context of hybrid 
schemes, it is customary to specialize the system (9)- (14) to the two- fluid case and to neglect the 
inertia of one of the fluid species (electrons). This last approximation is equivalent to taking the limit 
m2 for the second species in the total fluid momentum equation. Under this assumption, the 
sum of the equations (9) for s = 1, 2 produces 

pi-^ + Pi {Ui-V)Ui = (aipi + a2P2)'E + {aiPiUi + a2P2U2) x B - Vpi (23) 
Also, upon assuming neutrality by letting eo — >■ 0, the electromagnetic fields satisfy the equations 

Va,p,J7, = — VxB- /v/d^v, (24) 
^ Po J 

- = -VxE, (25) 

Y,(^sPs = -qh j f d\ , V-B = 0. (26) 
s 

Then, equation (23) becomes 

p^ + p{U ■V)U = -{qh jfd\^B+(^-^V xB- ah j-^fd\^ x B - Vp , (27) 

where we have dropped labels for convenience. Finally, inserting Ohm's ideal law E -|- [/ x B = 0, the 
kinetic two-fluid system becomes 

p^ + p{U ■V)U = {qhU j fd\- g,,yv/dV + -i-VxB^ x B - Vp (28) 

|^ + V-(p[/) = (29) 

^+v-|^ + a,(v-J7)xB.|^ = (30) 
ot ax av 

^ = V X ([/ X B) . (31) 

This is the same as the current-coupling hybrid scheme presented in [11, 36, 1], except that particle 
dynamics is governed by the Vlasov equation rather than its gyrokinetic counterpart. Notice that 
the above system does not make any assumption about the form of the Vlasov distribution for the 
energetic particles. Therefore, this system should in principle apply to a variety of other possible 
physical situations, as well. 

3.2 Euler-Poincare reduction by symmetry 

We now turn our attention to the Euler-Poincare formulation of these equations. That is, we ask 
whether the above current-coupling system possesses an Euler-Poincare variational principle. A pos- 
itive answer is provided by the reduced Lagrangian 

1{U, p, X, /, A) = i y p\U\' d^x - j pU{p) d^x - ^ y I V X Aj^ d^x 

+ y'/(^|u|2 + !|i|u-v|2 + g,(u-[/).A)d3xd3v (32) 
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-- 0, 


dA 

~dt 


d SI 
dt 


51 
5X 


= /V(: 




= 0, 





of the type 

together with the fohowing Euler-Poincare theorem. 

Theorem 2 The hybrid current- coupling MHD scheme (28) -(31) arises from the Euler-Poincare vari- 
ational principle 

■ ['l{U,p,X,f,A)dt = 

with the Lagrangian in (32) and variations given by 

5U = dtW-£uW, 5X = dtZ-£^Z, 6f = -£xf , Sps = -£uP, SA = -£uA 

where W and Z vanish at the endpoints. This variational principle is equivalent to the Euler-Poincare 
equations 

" " - - ^-4) A, (33, 

£uA = 0, (34) 

(36) 

which hold for an arbitrary hybrid Lagrangian. 

Proof. The derivation of the Euler-Poincare equations (33)- (36) from the Euler-Poincare variational 
principle can be easily obtained by direct verification [21, 7]. In order to derive the current-coupling 
MHD scheme (28)-(31), one simply computes the functional derivatives and inserts them into the 
Euler-Poincare equations (33)-(36). In particular, for the Vlasov kinetic part one has 

SI „, ^ 51 51 mh, ,9 mh , ,n 

— =mhfi2u-v + ahA), ^=0, — = u ^ + u - v ^ + g?,, u - [/ • A 

ou oa of Z z 

Then, on projecting equation (18) onto the second component (recall that X = (u, a)), we get 

= 5vU • d^— = u(x, v) - V 

/ Sn Of 

so that 

X(x,v) = (v,a(x,v)) , — = mhf{yr + ahA), — = ^jvp + g;, (v - (7) • A . 

ou of I 

Upon denoting p(x, v) = m/^v + q^ A(x) and dividing equation (18) by /, one finds 



which when projected onto the first component yields 

dtp + (v • 5x + a- av)p + (9xv - p = 9x(^^|vp + qh{-v - U) ■ a] . 
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Upon recalling that (x, v) are independent coordinates and using standard vector identities, we can 
write 

a(x,v) = -ahiV{U ■ A) + ^ j +afev x (V x A) 
= a,, (v - ?7) X (V X A) 

where the bottom line is justified by the second equation in (34). Therefore, equation (36) returns 
the Vlasov kinetic equation (30) in the form 

| + v.g + .4(v-t/)x(VxA)].g = 

with a modified Lorentz force. 

We now focus on the fluid part. It suffices to compute 

^=pU-qhnA, ^ = -V xV xA + qh{K-nU), = ^\uf + pU'{p) +Uip) 

where we have introduced the additional notation 

n = jfd^v, K = jvfd\. 

At this point, it suffices to insert the above functional derivatives into equation (33), so that 

+ "^t/j (pU -qhnA)=p\/ Q \U\^ + pU'{p)+Uip)^ + (V x V x A) x (V x A) 

-qh{K-nU) xV X A + QhAV ■ (K-nU) (37) 

We observe that the zero-th moment of the Vlasov equation (30) satisfies dtn + V • K = 0. Then, 
making use of the second equation in (34) yields 

{dt + £u){nA) = -AV ■ (K-nU) (38) 

while expanding the Lie derivatives in (37) returns the velocity equation 

p^+p{U ■V)U = (qhuU- qhK + —VxB) X B - Vp (39) 

in which we have substituted B = V x A and p = p^hl'{p). ■ 

Remark 4 (Euler-Poincare reduction) Upon following the treatment in [21, 7], one finds that the 
unreduced Euler-Poincare Lagrangian of the current- coupling scheme is a functional of the type 

Lpo,AoJo ■ TDiff(M3) X rDiff(M'5) ^ M 

where x denotes direct product. Consequently, 

The reduced Euler-Poincare Lagrangian (32) is obtained by the reduction process 
l{u, X, p, A, /) = Lpoo,,-!, Aoor,-i,/oov-i (?) o -0 o i;-^) . 

Here the advected quantities p, A, f are acted on by the corresponding diffeomorphism groups, taking 
into account their intrisinc tensorial nature, that is (p. A, /) G Den(R^) x f]^(]R'^) x Den(]R^), where 
n^(M^) denotes the space of differential one-forms on M^. 
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3.3 Discussion 

Kelvin circulation law. Relation (37) amounts to the following Kelvin circulation law 

— / (U -Qh-A] ■dx = qh J) - ({V ■ (K-nU))A- (K-nU) xb) -dx 
Jy,{u) ^ P ' Mu) P ^ ' 

+ (f - (V X B) X B • dx, (40) 

which agrees with the corresponding result found in [44]. Notice that the creation of circulation 
on the right hand side is generated by the terms involving 61/5 A in the Euler-Poincare equation 
(33). As explained in [21], these terms comprise a momentum map generated by the action of the 
diffeomorphisms on the cotangent bundle r*r2^(]R). The presence of these terms is related to the 
fact that the non-zero magnetic potential A (together with the mass density p) breaks the relabeling 
symmetry of the unreduced Lagrangian, so that 

On the other hand, the Kelvin circulation theorem for the hot particles reads simply 




p • dx = , 



which recovers the well known preservation of the Poincare invariant for the hot particle motion. 



Ertel's theorem. The above Kelvin circulation law identifies the expression of the force 

* = (V • (K - n J7)) A - (K - n ?7) X B + /i(7^(V X B) X B (41) 

acting on the fluid with momentum pU — QhnA. The above quantity can be used to generalize 
Ertel's theorem for MHD (see [15] and references therein) to the hybrid current-coupling scheme. For 
simplicity, consider the incompressible case, so that p = 1 enforces V • ?7 = 0. Next, project the 
quantity p~^nA onto its divergence-free part by defining [p~^nA] = p~^nA + Vip, for a scalar 
function if such that V • [p~^ n A] =0. Notice that we keep the density p in these relations to provide 
correct dimensions, while p = 1 for incompressible flows. Then, upon denoting Dt = dt + U ■ 'V and 
w = V X (U — Qh \p~^ n a]) it is easy to see that the curl of equation (37) produces the generalized 
Ertel relation 

Dt [Q ■ Va) - (cj • V) Dta = Vq • V x * , (42) 

where a is an arbitrary smooth function and ^ is the force expressed by (41). The quantity lj ■ Va 
is the potential vorticity and the above relation generalizes Ertel's theorem to the current-coupling 
scheme of hybrid MHD. 

Cross helicities. Upon denoting V = U — qhP~^ n-A, the following two cross-helicities may now be 
defined: 

Ai = / J7 • B d^x A2 = / V • B d^x 



However, while the first is conserved in time, i.e. dAi/dt = 0, the second satisfies 

^A2 = -g^^ j p-^ nA ■Bd^x = I p-' {A ■B)V ■ (K - nU)d^^, 
where the last non- vanishing integral is generated by the term parallel to A in equation (37). 
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4 First pressure-coupling hybrid MHD scheme 
4.1 Formulation of the model 

In this section we show how the variational structure of the previous current-coupUng scheme provides 
a basis for the Euler-Poincare formulation of a pressure-coupling scheme. This scheme establishes an 
equation for the total velocity 



U = U + — / v/dV, 
P J 

under the assumption that the kinetic moment 

K= [vfd^v (43) 



does not contribute to the total energy of the system. This assumption can be justified if the energetic 
component is particularly rarefied, so that its density 

n = y / d^v (44) 

is negligible compared to the density p of the cold fluid. This is precisely the hypothesis that we shall 
use in our derivation of the following energy-conserving pressure-coupling scheme obtained in [44] 

+ p(U ■ V)U = -Vp - m^V • / vv/d^v - — B x V x B (45) 
ot J po 

-wr + iU + v)-- — -VU -v + a^v xB ■ — =0 (46) 

ot ox Ov av 

^ + V.(pU) = 0, ^ = Vx(t/xB). (47) 

Here we have dropped the bar symbol for convenience. Before proceeding further, we remark that the 
L/"-terms appearing in the kinetic equation (46) differ substantially from the corresponding term in 
the hybrid MHD model presented in [30] (whose Vlasov kinetic equation is replaced by its gyrokinetic 
approximation). Indeed, the fluid transport term U ■ d^f is totally absent in reference [36], where the 
circulation force term d-^f ■ VU ■ v is replaced by the Lorentz force —QhU x B emerging in (30) as an 
electric field force (from ideal Ohm's law). More particularly, upon denoting IP = / vv / d'^v, reference 
[36] derives a pressure-coupling scheme by assuming S^K = —V • P -|- a/j(K — nU) x B ~ in the 
current-coupling model (28)-(31); then the resulting force balance allows replacing Lorentz forces by 
the pressure term in the momentum equation (28). These crucial steps break the energy-conserving 
nature of the system, as explained in [44]. However, notice that the static equilibria of the above 
equations (45)-(47) coincide with those of the hybrid model in [36], provided the hot particles are 
governed by Vlasov dynamics. 



4.2 Euler-Poincare reduction by symmetry 

Although the physical approximations leading to the pressure-coupling scheme present some problems 
that were summarized in [44], we shall see below how the variational approach to the model (45)-(47) 
produces an Euler-Poincare system on a semidirect-product Lie group. At the reduced level, we shall 
prove that the Euler-Poincare Lagrangian is a functional of the form 

/ : (je(]R^)(DX(]R^)) X Den(]R3) x n^(R^) x Den(M^) ^ R 



D.D. Holm & C. Tronci 



Hybrid plasma models 15 



where the infinitesimal action that is involved in the semidirect-product Lie algebra X(M^) (S)X(M^) 
is given by 

U-X = £xu^, where X[7 := (?7, (v • V)?7) V J7 G X(M^) . 

This action naturally arises from the tangent-lifted action of Diff (M'^) on TW^ = M^, which in turn 
generates the natural Diff (M^)-action on Diff (M^) (see Remark 5 below). On the other hand, the space 
of the advected quantities [21] 

(p, A,/) € Den(M2) x n'^(M.^) x Den(]R^) 
involves the Lie algebra representation 

((7,X) • (p,A,/) = {£up, £uA, £x+xuf) 
whose associated diamond operation, defined by 

{{rA-Jf)'' ^' ^' ■ - - ( (I' ^' I?) ■ (^^ ^' ^) > ■ 

will be derived explicitly in what follows. 

At this point, the problem has been cast into the standard Euler-Poincare theory for parameter- 
dependent Lagrangians L : TG x V* ^ M, with the peculiarity that the Lie group G is a semidirect- 
product. Indeed, upon replacing G by G(S)-ff, the pressure-coupling scheme will be written as an 
Euler-Poincare variational principle on T{G(§)H) x V*. Upon specializing to the case G = Diff(]R'^), 
H = Diff(r]R3) and V* = Den(R3) x 0^(]R3) x Den(R6), the Euler-Poincare equations associated to 
such a Lagrangian can be written as follows on the reduced space jC(]R^) (S)X(M^) x V*: 



9^ ^ + ^ = ^ * ^ + ' ^' ^ J ^' ^) 

^ + £up=0, ^ + £uA = (50) 
d SI ^ 61 61 

^^ + i:x+x,^ = /V(.,.)- (51) 
df 

+ £x+Xu f = (52) 

where (oi) in equation (49) denotes the Z7-component of the diamond operation defined in (48), and 
the star (*) operation is defined as 



5X ' / \5X 



★ X, J7 ):=-(—, i^XE/X). (53) 



Integration by parts yields the more explicit expression, 
61 \ I 61 

i^x^, (t/,(v-V)[/) 
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so that 



where the indices 1 and 2 denote the u- and the a-components, respectively. 

In order to complete the set of equations (49)-(52), we shall need a suitable Euler-Poincare La- 
grangian, which is given by 

1{U, p, X, /, A) = i y p\Uf d^x - I pU{p) d^x - ^ y I V X A|2 d^x 

+ |/(^|u|V^|u-v|Vg.u.A)d3xdV (55) 

This Lagrangian is obtained from (32) by simply neglecting the term q^JnU- Ad'^x, consistently 
with the assumption of a rarefied energetic component. 

The simplest starting point involves the kinetic part of the system (49)-(52), which is composed 
of the last two equations, i.e. (51)-(52). Let us start by calculating the functional derivatives. Upon 
using similar arguments as those in the previous section (and especially using the second component 
of equation (51)), in slightly different notation X(x, v) = (u(x, v), q:(x, v)) one finds 

X + Xt7 = (v + J7,a + (vV)t/), ii = „i,/(v + a,A), = !I^|v|2 + • A . 

' OM OJ 2 

Upon denoting p(x, v) = m^v + A(x), we divide equation (18) by / so that 

d /I 6l\ „ / 1 SI 



and project it onto the first component to obtain 

dtp+{{v + u)- d^)p + {cx ■ d^)p + ((v • d^)u ■ 5v)p + (v + J7) • p = ax(^|v|2 + ^^^v a 

Then, upon using the second equation in (50) as well as standard vector identities, one writes 

a = a/jV X (V X A) - VU • v - (v • V)f/ 

and the vector field X + Xfy becomes 

X + Xj7 = (v + i7, ah^ X (V X A) - VJ7 • v) . 

In turn, upon noticing that V(x,v) " (^ + Xjy) = 0, this vector field produces the Vlasov kinetic 
equation (51) in the form 

§[ + (v + [/).|^-(V[/.v-a,vx(VxA)).|^ = 

which is identical to (46). 

At this point, one needs to verify that equation (49) effectively returns the velocity equation (45) 
of the pressure-coupling scheme (45)-(47). To this purpose, we shall use the following 



Lemma 3 In the special case when 

fx. vn n.n.d 

6X 



51 

X(x, v) = (v,q:(x, v)) and — (x, v) = (w(x, v), O) 



then 

for arbitrary vector quantities Q;(x, v) and w(x, v). 
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Proof. The proof follows by direct verification, upon writing the definition of the star operation in 
(53) as follows 



(5X' / \(^X' 



51 



— , y (X}7 • V(x,v)) X - (X • V(x,v)) Xj7 

^w, (Xj7 • V(x,v)) V - (X • V(x,v)) U 
|w, (((v • V)U ■ 5v)v - (v V) (7 

|w, (^(v • v)(7 - (v • V) jy) ^ = 

in which the last step uses integration by parts of the first term. ■ 

Another result that we shall need is the following formula for the diamond operation in (53) 

fjOif = jfd^fjd^-V.jf^rd^fjd^, (56) 
which may be directly verified from its definition 



61 \ / / -.61 

jjo, /, C/^ :=^/, (X,;.V(.,v))- 



Now, upon recalling the particular form of the variational derivative 

= -jjri'iA'^] + • A, 

we calculate 

1^ oi / = qhVA ■ K - rrihV • P - g/. (V • K) A - g;, (K • V) A 

= (7hK X B - (7;, (V • K) A - m^V ■ P , 

where we recall the definiton (43) of the averaged kinetic momentum and we have introduced the 
absolute pressure tensor 

P = jvvfd\. 

Then, upon writing 
and evaluating 

Jlj=P^^ ^ = -VxVxA + qnK, ^ = i \U\^ + pU'ip) + U{p) 

we see that equation (49) returns the velocity equation (45) of the pressure-coupling scheme (45)- (47). 
In conclusion, we have proven the following theorem. 

Theorem 4 The hybrid pressure- coupling MHD scheme (45) -(47) arises from the Euler-Poincare 
variational principle 



5 ri{U,p,XJ,A)dt = 
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with the reduced Lagrangian 

I : (X(M^)©X(M^)) X C°°(]R3)* X n'^{R^) x C°°(R^)* 
given in (55) and variations 

d{U, X) = dt{W, Z) - {£uW, £xwX -£^u'Z + £^Z) 

(^Z = --^z+Xw/' 6{p,A) = -£w{p,A.) 

where the vector fields W € X(M'^) and Z € X(M^) vanish at the endpoints. This variational principle 
is equivalent to the Euler-Poincare equations (49)-(52), which hold for an arbitrary hybrid Lagrangian. 



Remark 5 (Conjugation action in semidirect-product Lie groups) The Lie algebra action that 
is involved in the semidirect product X(M^)(DX(M^) is naturally inherited from the Jacobi-Lie bracket 
on X(M^). According to the theory of semidirect-product Lie groups, this action must arise from a 
group action o/Diff(M^) on Diff(R^) that is also a group homomorphism. In other words, rj {ipi V'2) = 
r] (ipi) r] (ip2) , with rj € Diff(]R^) and V'i;V'2 G Diff(R^). In particular, since we can regard Diff(M'^) 
as a subgroup o/Diff(M^), one is led to consider the action ip 1-^ rj o ip o r]~^ which is naturally in- 
herited from the conjugation action in Diff(M^). This action generates the semidirect-product Lie 
group Diff(M^)(s) Diff(M^), whose tangent space at the identity (§)X(R^) is endowed with the 

Lie bracket 

[{U, X), {W, Z)] = -{£uW, i^XwX -£^u^ + £^Z) . 

More details on semidirect-products of two Lie groups may be found in [28, 4]- 

4.3 Discussion 

The Euler-Poincare construction of the first pressure couphng scheme is based on the following propo- 
sition: 

Proposition 5 The Euler-Poincare equations (45) -(47) yield 

Proof. The proof is a direct verification, based on relations (54) and (56). After computing 
/ B ^ '\ / SI fSl., [, „ ^ 51 \ SI «! . 

-/(v.aj(^x„^)^dV + ^,/(v.aj£dV, 

the proof follows immediately from Lemma 6 below. ■ 

Kelvin-Noether theorem and its momentum map. The above relation represents the La- 
grangian analogue of an important construction in Lie-Poisson Hamiltonian systems, known as un- 
tangling. Untangling is accomplished by shifting the momentum by a momentum map that takes the 
Lie-Poisson bracket on the dual of a semidirect-product Lie algebra into the Lie-Poisson bracket dual 
to a direct-sum Lie algebra. For more details, see Corollary 2.4 in [24]. 
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It is perhaps not surprising that the very first application of this construction occurred in plasma 
physics [18, 17]. This construction was also used in [44]. The momentum map in the present case is 
the dual i* : X*(M^) X*(M^) of the Lie algebra inclusion i : U ^ X[/. The result (57) hinges on 
the following property, which is proven in Appendix A: 

Lemma 6 The following map i* : X*(]R^) 
is a momentum map satisfying the relation 




for an arbitrary vector field U € X(M'^). 

Notice that the momentum map i* : X*(]R^) — )> X*(M'^) is different in nature from the star operator 
-k : jC(]R^) X j£*(]R^) — > j£*(ffi^) introduced in (53). Indeed, while the latter arises from the cotangent 
lift of the Diff (M^)-action on 3£(M^), the momentum map i* arises from the Diff (M^)-action on Diff (M^), 
which is given by conjugation, as explained in Remark 5. The momentum map property of i* can be 
easily verified since the inclusion i : ^ X(M^) is the dual of a Lie algebra homomorphism, i.e. 

[Xt7,Xw] = m),iiW)] = i{[U,W]) = X[u,w] , 

where [•, ■] denotes minus the Jacobi-Lie bracket on X(M'^) or jC(]R^), depending on the context. The 
proof in Appendix A shows explicitly that i* satisfies the definition of momentum map. 

An immediate consequence of Proposition 5 is the following circulation law for the hybrid scheme 
(45)-(47), which recovers the previous results in [44]. 



Corollary 7 (Kelvin circulation law) The pressure- coupling MHD scheme (45) -(47) possesses the 
following equivalent circulation theorems 

— (f U ■dx = - (f - i—B X V X B + mhV ■ [ vv/ d^v ) • dx (58) 
J-yt{U) ht{u) P \Po J ) 

±1 fj7-i //pdO- dx= - / i(B X (^-iVxB-(7,,K) -(7,,(V-K)a) -dx. (59) 

Proof. Upon considering the Euler-Poincare Lagrangian (55), relation (58) is implied by the Euler- 
Poincare theorem. See [21]) and the equation of motion (49). On the other hand, relation (59) is an 
immediate consequence of equation (57) in Proposition 5. ■ 

Notice that taking the difference of the above circulation laws yields 

where we recall the relation p = m^v + A. Thus, upon considering (51) and the zero-th moment 
equation dtu + V • (nil) = —V • K associated to (46), we have 

/ — • dx = - / - f V • /"vv/ dV - UhK X bV dx . 

J-ft{U) P JytiU) P^ J ' 
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Moreover, the above fluid circulation laws are accompanied by preservation of the Poincare- invariant: 

d 



p • dx = , 

where the curve Ct now moves along the total phase-space vector field X + X[/ . 



Ertel's theorem. By proceeding as in the corresponding treatment for the current-coupling scheme, 
one recognizes that the force 

* = (//o'VxB-g;,K) xB + (7/,(V-K)A (60) 

provides the opportunity to generalize Ertel's relation for MHD [J 5] to apply to the first pressure- 
coupling scheme. Indeed, upon following similar steps as those in Section 3.3, one finds that the 
incompressible form of equation (57) yields the relation (42), with a) = V x (U — p~^J /pd'^v) and 
^ as given in (60). Again, we kept the density p in the expression of u) in order to provide correct 
dimensions; incompressible flows are always accompanied by p = 1. 

Cross helicities. Notice that, upon denoting W = U — p~^ f f pd^'v, both of the cross helicities 

Ai = j U -Bd^x A3 = W • B d^x 

possess nontrivial dynamics. Indeed, their equations of motion read as 

'^ = -mnj p-' (V • P) • Bd^x, ^ = (A ■ B) V • Kd^x. 

On the other hand, the following cross helicity is conserved: 

A2 = Jr-Bd^x, (61) 

where we have denoted T = U — m/i/9~^K. Upon noticing that Dt{p~^ n) = p~^{dtn + V • (nil)), the 
conservation of A2 is readily seen by computing 



dA2 dAa 



dt dt 



+ p~^nA-Bd^x = 0, 



where one considers the equation dtn + V • (nil) = — V • K arising from the zeroth moment of the 
Vlasov equation (46). 

5 Second pressure-coupling hybrid MHD scheme 
5.1 Formulation of the model 

As mentioned in the previous section, the pressure-coupling MHD scheme is conventionally obtained 
under the assumption that the hot plasma component is rarefied. Upon denoting by n the particle 
density of the hot component and by mc the cold particle mass, this assumption reads as n <C p/ruc- 
Then, in order to avoid divergences in the mean velocity V = n^^ f d'^v of the hot component, a 
small hot particle density n requires the hot momentum K. = fv f d'^v to also be small. Thus, it is 
customary to replace the total momentum pU + m/jK by simply pU, i.e. the cold fluid momentum. 
While this operation is often performed on the equations of motion [8, 36], our approach makes this 
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replacement directly in the variational principle, resulting in agreement with [44]. The advantage 
of modelling in the Lagrangian of the Euler-Poincare variational principle is that it always produces 
circulation theorems. This is the content of the Kelvin-Noether theorem of [21]. 

The assumption of a rarefied hot component may also require that the mean kinetic energy 
m/i/2 Jn|V|^d^x of the hot component is subtracted from the corresponding total kinetic energy 
mfi/2 J f jv|^ d'^xd'^v. This operation yields the second pressure-coupling scheme [44] 

+ p{U ■V)U = -Vp-mhV ■ j (v-^^ ^v- /dV- -^B X V X B (62) 

I i - - ^) ^ ■ - (v^) ■ I - ^)) ■ 1^ ^ » <») 

^ + V-{pU) = 0, ^ = Vx([/xB). (64) 

We remark that neglecting all U- and n~^K-terms in the kinetic equation (63) and replacing ?i~^K x B 
by ?7 x B produces the hybrid MHD model in [2.3, 43] (although the general Vlasov equation is adopted 
here, rather than a drift-kinetic equation). Upon denoting F = /(v — (v))® ^/d^v and (v) = K/n, the 
model in [23, 43] can be derived by assuming n 9t(v) + n ((v) • V) (v) = —V • P-l- a/j(K — nil) x B ~ 
in the current-coupling scheme (28)-(31), so that Lorentz forces in (28) are replaced by a relative 
pressure term. 

Notice that the static equilibria of the above equations (62)-(64) coincide with those of the hybrid 
model in [23, 43] (for hot particles undergoing Vlasov dynamics), provided the equilibrium Vlasov 
distribution (usually denoted by /o) is isotropic in the velocity coordinate, i.e. Kq = / v/od'^v = 0. 



5.2 Euler-Poincare reduction by symmetry 

In order to obtain the Euler-Poincare formulation of the hybrid model (62)-(64) for the second pressure- 
coupling scheme [44], the Lagrangian (55) is transformed into 

/([/, V, X, p. A, /) = i \Uf d^x - I pU{p) d^x - ^ y I V X d^x 

+ //(^^|u| + — |u-v| - — \V\ +g/,(u + F) • Ajd^xd^v, (65) 

where we notice that the mean velocity Y appears as a new dynamical variable. The term Qh J fV ■ 
Ad^xd^v has been inserted in order to match the correct Lorentz force on the hot component [4 1]. As 
we shall see, the two V— terms in the above Lagrangian correspond to subtracting the contributions of 
the mean velocity (v) = n~^K = —V. Upon following the same reasoning as in the previous section, 
we realize that the fluid [/-transport exerted by the cold fluid component on the Vlasov distribution 
/ of the hot particles must imply a [/-transport of the mean hot velocity V. More particularly, we 
interpret the above Lagrangian as a functional of the type 

I : jei(M^)(D(X2(M^)(DX(M^)) X Den(M^) x n\R^) x Den(M6) ^ M, 

where Xi(M^) and ^£2(1^^) are two copies of the same Lie algebra jC(M^) of vector fields, although they 
are denoted differently because the second is assumed to act trivially on the space Den(M3) x ^^^(IR^), 
containing the cold fluid density p as well as the magnetic potential A. The first (outer) semidirect- 
product symbol corresponds to fluid [/-transport of both the mean velocity V G X2(M^) and the 
phase-space vector field X € X(M^). On the other hand, the second (inner) semidirect-product symbol 
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corresponds to the V-transport exerted by the mean flow of the hot component on its corresponding 
phase-space velocity. At the group level, the unreduced Lagrangian is of the type 

Lpo,AoJo ■■ T(Difri(R3)(D(Diff2(M=')(S)Difr(M6))) ^R, (66) 

where {po, Aq, fo) are the advected parameters. Notice that similar arguments to those in Remark 5 
also apply here about the group actions involved in nested semidirect-product Lie group structures of 
this kind. The first instance of nested semidirect-product Lie-group structures also occurred in plasma 
physics: in the discovery of the Lie-Poisson brackets dual to nested semidirect-product Lie algebras 
in models of Alfven wave turbulence [14, 16]. This construction was also used for hybrid Vlasov-fluid 
models in [ , ,]. Further details can be found in [12], where similar Lie group structures were shown 
to arise in polymer dynamics. 

At this point, general geometric mechanics arguments ensure that the Euler-Poincare variational 
principle 6 f^^ 1{U, V, X, p, A, f)dt = produces the following equations of motion: 



d 61 61 61 SI ^ f6l 61 61 , . ^ ^. 

+ £c/ 777 = --^VTTT + TTF * X + — , — , — oi (p, A, /) (67) 



dt6U ^ 6U 6V 6y^ \6p 6A' 6f 

^ + £up=0, ^ + £uA = (68) 
d 61 ^ 61 61 ^ 61 ^ 

d 61 ^ 61 61 

+ i?x+x.+, ^ = / V(x,v)^ (70) 

^ + £x+x^^u f = (71) 

In order to see, how equations (67)-(71) recover the second pressure coupling scheme (62)-(64), it 
suffices to substitute the Lagrangian (65). After computing 

— = m/j/ 2u - v + Oh A) , ^ = 0, — = — u + — u-v - — \V\ +qh{n + V) ■ A, 
ou oa oj 2 2 2 

the second component of equation (52) yields u = v so that X = (v, Q;(x, v)), similarly to the results 
in the previous section. Moreover, the first component of (52) reads as 

9tp + ( (v + y + i7) • c>x)p + (a • d^)p + 5x(v + y + (7) • p = ^^((//.(v + y) . A - ^jFp 

where we have denoted by a = ct + (X.y+u)2 = ct + (v • V)("\^ + U) the total acceleration of the hot 
particles. After using standard vector identities, this equation yields the expression for the total force 
on the hot component: 

m/ia = Qh + V) X B - ruhVU ■ v — mhW ■ {v + V) . 

Then, the total vector field X + 'Kv+u = {v + V + U,a) is divergence-less and the Vlasov equation 
(52) becomes 

^ + (v + V + (7)-|^ + (a;,(v + V)xB-V?7-v-VV-(v + F)) -1^ = 0. (72) 

At this point, we observe that the constraint V = — K/n is preserved by the dynamics. This is a 
direct consequence of the following 
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Proposition 8 With the notation of Lemma 6, the Euler-Poincare equations (67)-(71) yield 



dt ^^"^ J \5V \6X^ , 

The proof proceeds analogously to that of Proposition 5 (see also Lemma 6). Then, upon considering 
the Lagrangian (65), the Euler-Poincare equations (67)-(71) preserve the constraint 

-V jfd^v = jvfd^v, 

which allows one to recover the Vlasov equation (63) of the second pressure coupling scheme in [44]. 
Analogously, one can show that equation (67) recovers the hybrid equation of motion of the same 
pressure coupling scheme, that is (62). In order to show this, it suffices to verify that 



This formula requires a lengthy but straightforward calculation that uses V = — n^^K and the well 
known relation 

V-y"vv/dV = V - (^n^^KK + j (v - n'^K) (v-n-^K)/dV^ 
between the absolute and relative pressure tensors. Thus, in conclusion, we have proven the following. 



Theorem 9 The hybrid pressure- coupling MHD scheme (62) -(64) arises from the Euler-Poincare 
variational principle 

6 [\{U,V,X,p,AJ)dt = 

J to 

with the reduced Lagrangian 

I : Xi(M3)(s)(X2(M^)(S)X(M6)) X C°°(M3)* x n\M.^) x (:7°°(M6)* 

as in (65) and variations 

6{U, V, X) = dt{W, P, Z) - {£uW, £wV - £uP + £v'P, i^Xp+w^ - i^x^+^Z + £xZ) 
^/ = -"^^z+Xp+w/ ' '5(p, A) = -i^w (p, A) , 



in which the vector fields P,W € X(]R ) and Z G X(M ) vanish at the endpoints. This variational 
principle is equivalent to the Euler-Poincare equations (67)-(71) which hold for an arbitrary hybrid 
Lagrangian. 

5.3 Discussion 

It is relevant to notice that equations (67) and (70) yield the following relation: 

61 / 5l\ 

__xVxA+(v.^)a^ (74) 



D.D. Holm & C. Tronci 



Hybrid plasma models 24 



where one has used Lemma 6. Upon inserting the Lagrangian (65), the Lie derivatives in the right 
hand side cancel due to the constraint V = — K/n. 

Kelvin circulation laws. The Kelvin-Noether conservation laws also hold for the equations (62)- 
(64). Indeed, upon repeating the same steps as in the proof of Corollary 7, one finds the circulation 
laws 

/ U-dx = -(f i ( — B X V X B + m/iV • P ) -dx (75) 
J-ytiU) J-yt[U) P \fJ-o ) 

A/ ( ?7- i //pd^v ) •dx = -/in^ / -BxVxB-dx, (76) 
J^t{V) \ P J ) htiU) P 

where the second is a direct consequence of equation (74). These results coincide with those found in 
[44] within the Lie-Poisson Hamiltonian setting. Now taking the difference of the above two relations 
yields 

it/ ( - /"/P d^v ) • dx = m/, ^ / — • dx = -m/, / i(V-P)-dx, 

di J-ttiU) \PJ J dt /^^(t/) p J^,(u) P 

where we have used the zero-th moment equation dtn + V • [nU) = associated to (72). Indeed, 
together with equations (68), this yields the following conserved circulation: 

A/ + A.dx = 0. 

di/7i(j7)V PJ 

Moreover, the above circulation laws are accompanied by the following Poincare invariant relation: 

4 / p • dx = , 

dt J<t(x+Xcj_K/„) 

where the curve Qt moves along the total phase-space vector field X + Xjy_K;/„. 

Ertel's relation. Proceeding as in Section 3.3, taking the curl of the incompressible version of 
equation (74) for the Lagrangian (55), or equivalently applying the Stokes theorem in the Kelvin- 
Noether relation (76) produces an Ertel relation of the form 

Dt (cj • Va) - (cj • V) Dta = -^q ^ Va • V x (^B x (V x B)) , 

where a is an arbitrary scalar function and, upon keeping the mass density p (equal to one for 
incompressible flows), 

= V X \ U — mh — 
V P 

The Ertel relation above written in terms of u> has the same form as the corresponding relation in 
[15], except that the vorticity Q here involves a velocity shift due to the mean specific momentum 
carried by the particles. 

Cross helicities. Upon denoting W = U — p~^ f f pd^v, it easy to see that the two cross helicities 

Ai = [ U -Bd^x and A3 = /" W • B d^x 
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possess the following dynamics 



dt J ^ ' dt 

so that A3 is now conserved by the hybrid dynamics of equations (62)-(64). Also, notice the conser- 
vation of the following modified magnetic helicity 

n= /-A-Bd^x, 

J p 

which allows to write the cross helicity invariant A3 as 

A3 = y" {u-nih^^ -Bd^x 

The conservation law for % provides an interesting opportunity to study the stability properties of 
this hybrid scheme. In particular, % does not vanish for static equilibria. This means the energy- 
Casimir method may be applied for hybrid fluid equilibria that are analogous to the Chandrasekhar 
flows of inviscid MHD, [221. 



6 Summary and conclusions 

This paper has derived three different hybrid Vlasov-fluid plasma models by using the Euler-Poincare 
approach first developed for the Maxwell- Vlasov plasma [7] . After presenting the Euler-Poincare ap- 
proach for the Vlasov-multifluid plasma system, the discussion focused on three different schemes for 
deriving hybrid Vlasov-fluid MHD models. These comprised the current-coupling scheme and two 
pressure-coupling schemes. The first hybrid model was written on the direct product of two different 
diffeomorphism groups, as explained in Remark 4. The second one involved the more sophisticated 
construction of the semidirect-product diffeomorphism group discussed in Remark 5. Finally, a com- 
pound semidirect-product structure arose for the third hybrid model, see equation (66). In all three 
theories, Kelvin circulation theorems were presented for both the fluid motion and the hot particle 
dynamics on phase space, and the invariant cross-helicities were identified. A Legendre transform in 
each case would recover the Lie-Poisson results found [11]. Shifting to the drift-kinetic (or even gy- 
rokinetic) approximation would require another Lagrangian, which may also be derived systematically 
from the Lagrangian for Vlasov- MHD. Summaries of the properties found here for the three different 
hybrid Vlasov-fluid plasma models are given below. 



Current- coupling MHD scheme 
• Equations of motion (28)- (31): 



+ p{U ■V)U = (qhuU- g,,K+— VxB) x B - Vp 
at \ fio 

^ + v.|i+«.(v-C/)xB.|i=0 

dt ax av 
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Kelvin circulation laws derived from (38)- (40): 

/ U ■dx= J) - ( QhuU - qhK + —V x B 1 x B • dx , 

at J^t Jjt P \ Po 

^ I (l + - ) A-dx = / i(v-(nJ7-K))A-dx 



<^tht\ PJ htP 
where 7t is any closed loop that moves with the fluid velocity U . 

• Magnetic and cross helicity invariants: 

?^ = / A • B d^x , A = / J7 • B d^x 



Approximation: same as ideal MHD; this yields the Lagrangian (32) 



Pressure- coupling MHD scheme - first variant 



Equations of motion (45)-(47): 



+ p(U ■ V)U = -Vp - m/jV • P - — B X V X B 
at fio 

| + (f/ + v).§^-|^.Vf/.v + a,vxB.|^ = 
ot ax ov ov 

_^ + V-(p[/) = 0, _ = Vx([/xB). 



Kelvin circulation laws (58)-(59): 



^ i U ■dx = - ([ - ( — B X V X B + mhV • P ) • dx 

dt y^,, P \Po 

— / — • dx = / i (uhK X B - V • P ) • dx 

dt J^^ p J^^ p 

d_ 
dt 



/ ( 1 + - ) A • dx = - / i (V • K) A • dx ; 

J-rt V Pj At P 



'7t \ ■J ft 

Magnetic and cross helicity invariants: 



?^=/A-Bd^x, A= / ( J7-m;,— ) •Bd'^x 



P 

Approximation: neglects minimal coupling term JnU ■ Ad'^x in the Lagrangian (32). 
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Pressure- coupling MHD scheme - second variant 



Equations of motion (62)-(64): 

+ p(U ■ V)U = -Vp -mhV -F- — B x V x B 
dt po 

§^ + (v+c/-^i).|^ + L(v-^^)xB-v^.v+(v^i).(v-^i)).|^ = o 

at \ n J ox \ \ n J \ ^ J \ n JJ ov 

do (?B 

J^ + V-{pU) = 0, _ = Vx(J7xB). 

Kelvin circulation laws (75)-(76): 

— / U ■dx = - (f - ( — B X V X B + mhV • P ) • dx 

^ / ^S.dx = - / i(V-P) - dx, ^ / + A-dx = 

Jy, P J-yt P dt J^^ V PJ 

Magnetic and cross helicity invariants: 

7^ = /( 1 + - ) A-Bd^x, A= /"(?7-m/,— ) -Bd^x 



PJ J \ P 

Approximation: neglects JnU ■ Ad^x as well as mean flow terms in the Lagrangian (32) 



The Euler-Poincare approach provided the means of comparing the geometrical properties of these 
three hybrid Vlasov-fluid plasma schemes in the same framework. This framework allowed the iden- 
tification and comparison of the geometric relationships within each scheme that were shared by the 
others. We expect that this framework will be useful in other modelling contexts. For example, one 
may imagine using the Euler-Poincare framework, (i) in the comparison and selection of Vlasov-fluid 
hybrid models, (ii) in the validation of previous derivations, (iii) in making choices among the schemes 
in various physical regimes, and (iv) as a basis for performing other derivations obtained by modelling 
in the Lagrangian. From the physical viewpoint the roles of heat exchange and other advected quan- 
tities should also lead to interesting effects in future investigations. For example, the introduction of 
another advected quantity would produce an explicit Ertel theorem for the evolution of potential vor- 
ticity. Of course, Ertel's theorem is an immediate result of the Euler-Poincare theory for any hybrid 
fluid- Vlasov model with advected quantities. However, it was discussed here only in the simple case 
of incompressible hybrid fluid flows. 

In another direction for future research in the context of potential vorticity, one may use Ertel's 
theorem to investigate the evolution of the gradient of the potential vorticity, as studied recently in 
geophysical fluid dynamics in terms of the vector B = VQ{q) x Va, [13]. Here Q is an arbitrary 
function and q = u: ■ Va is the potential vorticity arising in a given fluid model. Upon considering an 
advected function a such that dta + U ■ Va = 0, Ertel's theorem was shown in [13] to produce the 
dynamics of the vector 3 in the form 

dB 

dt ^ 
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where U = U — (aV x ^) and V x $ = Va x V (g Q' V • W) . Evidently, the divergenceless vector 
X> = V X $ breaks the frozen-in condition of the vector B and it thus affects the stretching properties 
that are governed by the left-hand side of the equation for the vector B, see [13]. The identification of 
the vorticity a) and the force ^ in each of the hybrid models discussed here would provide interesting 
opportunities to study the dynamics of the gradients of the potential vorticity. This would perhaps 
lead to the production of fronts and other fine structures at high wave numbers that develop by 
stretching of the B- vector. 
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A Proof of Lemma 6 



The map i* is easily seen to be a momentum map arising from the dual of the Lie algebra inclusion 
i : ?7 I— > Xj7. Upon denoting H = (5//(5X, the momentum map property 

{F,(i*(H),[/)}^ j,e)=t/.F(X*{M6))[i^] 



{F,(i*(3),[/)}. 



can be verified explicitly. Here, f/jr^^.^ige)) [F\ denotes evaluation on the functional -F(S) of the 
infinitesimal generator t^j:-{x*{K6)) of the Diff(R^)— action on the space of functionals J-'(X*(M^)) on 
the one-form densities in jt*(]R^). Upon using the right Lie-Poisson bracket on j£*(M^), one computes 



£ 



= Ujr(^x*(me'))[F] ■ 

where [•, •] denotes minus the Jacobi-Lie bracket on X 
The rest of the proof proceeds in two steps. First, 

£^u^ d3v = |(i7-9x + (v-V)(7-5v)^d3v+ / ( V(7 • ^ + ax((v • V)J7) • 



6l_ 
5cx 



+ (V-t/ + 5v-((v-V)C/) )^)dV 



: !{u-d^)^ d\ + /((v v)u ■ av) dV + / ((v • v)vu) ■ ^ dV 

J on J du J ' da 

+ /(^'^ ■ I; + <^ I;) ''''' " J ■ ■ I;) ''''' 
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One also has 

/ (v • a.) (i:x, = |(v • 5.) . a. + (v . v)[/ . av) d^v 

+ y"(v • ax) (av((v • v)u) • dV 

+ /(v • Sx) (^(V • [/ + av • ((v . V)(7) ) £^ dV 

Then, for each term, one computes 

J (v . ax) (([/ • ax) ^^d'v = u■v|{^r. ax)^ dV + y Tv((vc/)(vax)) £ dV , 
|(v • ax) ((av((v • v)c/)) d^v = v[/.|(v • ax)£ d^v + |((v • v)vj7) • £ d^v 
/(V . ax) ((V . u) ^) d3v ^ / ((V . V) (V . ^) ) £ d3v + (V . ^) /( V. ax) £ d3v 

and 

y (v ■ ax) ((a, • ((v . v)?7)) d^v = ax • y"(av • ((v • v)?7))v£ d^v 
= -v- y ((vv)c/)£d3v-y(vax) (((^•v)c/-av)£) d^v 

= - /((v • V)V[/) . £ d^v - I Tr ((V[/)(vax)) £ d^v - ^(v • ax) (((v • V)[/ • av) d^v . 
Thus, in conclusion 

which completes the proof. ■ 
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